function br = isParamsGood_p(pf, pStr2Double, Opt)
    global Opt; 
    
    t0 = 10000;
    ting = 365;                
    
    resParams0 = CalcResParams(pf, pStr2Double, 0);                   
    sol = NNsysExact_ing0_p(t0, resParams0);    
    v1 = deval(sol, t0);
    vopt = deval(Opt.optsol0, t0);
    v1total = RNA2RNATotal(v1(1), v1(10), v1(9));
    vopttotal = RNA2RNATotal(vopt(1), vopt(10), vopt(9));
    
    d0 = abs(v1total - vopttotal);
    dcomp = vopttotal * 0.3333;

    resParams15 = CalcResParams(pf, pStr2Double, 1.5e+6);        
    sol = NNsysExact_ingNonzero_p(ting, resParams15, vopt);    
    d15 = deviationlog10(sol, Opt.optsol15, vopt, Opt.statConc, ting);

    resParams75 = CalcResParams(pf, pStr2Double, 7.5e+6);
    sol = NNsysExact_ingNonzero_p(ting, resParams75, vopt);    
    d75 = deviationlog10(sol, Opt.optsol75, vopt, Opt.statConc, ting);
    
    resParams150 = CalcResParams(pf, pStr2Double, 15e+6);
    sol = NNsysExact_ingNonzero_p(ting, resParams150, vopt);    
    d150 = deviationlog10(sol, Opt.optsol150, vopt, Opt.statConc, ting);

    br = [ (d0 < dcomp) (d15 < 0.27) (d75 < 0.27) (d150 < 0.27) ];  
end